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Abstract 

We  study  equilibria  of  the  Smoluchowski  equation  for  rigid,  dipolar  rod 
ensembles  where  the  intermolecular  potential  couples  the  dipole-dipole 
interaction  and  the  Maier-Saupe  interaction.  We  thereby  extend  previous 
analytical  results  for  the  decoupled  case  of  the  dipolar  potential  only  (Fatkullin 
and  Slastikov  2005  Nonlinearity  18  2565-80;  Ji  et  al  Phys.  Fluids  at  press; 
Wang  et  al  2005  Commun.  Math.  Sci.  3  605-20)  or  the  Maier-Saupe 
potential  only  (Constantin  e/ fl/ 2004  Arc/i.  Ration.  Mech.  Anal.  174  365-84; 
Constantin  et  al  2004  Discrete  Contin.  Dyn.  Syst.  11  101-12;  Constantin 
and  Vukadinovic  2005  Nonlinearity  18  441-3;  Constantin  2005  Commun. 
Math.  Sci.  3  531-44;  Fatkullin  and  Slastikov  2005  Commun.  Math.  Sci. 
3  21-6;  Liu  et  al  2005  Commun.  Math.  Sci.  3  201-18;  Luo  et  al  2005 
Nonlinearity  18  379-89;  Zhou  et  al  2005  Nonlinearity  18  2815-25;  Zhou  and 
Wang  Commun.  Math.  Sci.  at  press),  and  prove  certain  numerical  observations 
for  equilibria  of  coupled  potentials  (hetal  Phys.  Fluids  at  press).  We  hrst  derive 
stability  conditions,  on  the  magnitude  of  the  polarity  vector  (the  hrst  moment 
of  the  orientational  probability  distribution  function)  and  on  the  direction  of  the 
polarity.  We  then  prove  that  all  stable  equilibria  of  rigid,  dipolar  rod  dispersions 
are  either  isotropic  or  prolate  uniaxial.  In  particular,  all  stable  anisotropic 
equilibrium  distributions  admit  the  following  remarkable  symmetry:  the  peak 
axis  of  orientation  is  aligned  with  both  the  polarity  vector  (hrst  moment) 
and  the  distinguished  director  of  the  uniaxial  second  moment  tensor.  The 
stability  is  essential  in  establishing  the  axisymmetry.  To  demonstrate  that  the 
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stability  is  indeed  required,  we  show  that  there  exist  unstable  non-axisymmetric 
equilibria. 

Mathematics  Subject  Classification:  35Kxx,  70Kxx 


1.  Introduction 

The  Smoluchowski  equation  describing  the  behaviour  of  rigid  rod  ensembles  in  a 
viscous  solvent,  often  called  rigid  nematic  polymers  or  rod-like  dispersions,  has  recently 
received  significant  attention  from  analytical  perspectives  [4-7,9-13, 15, 16, 18, 19].  The 
Smoluchowski  equation  provides  more  accurate  descriptions  based  on  physical  principles 
than  the  mesoscopic  orientation  tensor  equation  with  approximation  closure  rules.  It  is 
mathematically  much  more  challenging.  Led  by  Constantin  et  al  [4-7],  the  mathematical 
groups  started  to  study  the  Smoluchowki  equation  directly.  Whereas  most  treatments  have 
investigated  equilibria  with  only  the  Maier-Saupe  interaction,  a  couple  of  papers  deal  with 
rigid  rod  ensembles  subject  to  only  the  dipole-dipole  interaction  without  the  Maier-Saupe 
interaction  [11,18].  Based  on  the  formalism  of  Doi  and  Hess,  a  new  kinetic  theory  for  imposed 
kinematics  of  magnetic  rod  suspensions  was  recently  proposed,  and  numerical  solutions  were 
calculated  in  simple  flows  [1,2]. 

Prior  to  the  coupling  of  an  external  hydrodynamic  or  magnetic  field,  it  is  prudent  to 
first  understand  the  space  of  quiescent  equilibria  of  the  Smoluchowski  equation.  When  the 
individual  rods  are  dipolar,  and  they  are  sufficiently  crowded  (semi-dilute)  as  to  be  affected  by 
the  Maier-Saupe  interaction,  this  means  that  one  must  combine  the  analyses  referenced  above 
and  address  the  coupling  of  the  dipole-dipole  interaction  and  the  Maier-Saupe  interaction. 
From  this  rigorous  basis,  one  can  then  entertain  imposed  external  fields.  It  is  noteworthy  that 
rigorous  proofs  for  the  Smoluchowski  equation  and  its  equilibrium  distributions  have  only 
been  provided  in  the  last  couple  of  years,  for  both  limits  of  decoupled  potentials,  even  though 
the  essential  results  were  known  by  formal  arguments  since  Onsager’s  seminal  paper  [17], 
and  modern  numerical  simulations  have  detailed  the  probability  distributions  for  quite  general 
interaction  potentials. 

This  paper  addresses  equilibria  of  the  Smoluchowski  equation  for  fully  coupled  potentials. 
We  are  guided  by  recent  numerical  investigations  [14],  for  which  we  prove  the  primary  physical 
results  regarding  the  properties  of  stable  equilibria.  In  particular,  it  was  shown  in  [14]  that  the 
polarity  vector  (first  moment  of  the  orientational  probability  distribution  function  (PDF))  of 
all  anisotropic  equilibria  must  be  a  principal  axis  of  the  second  moment  tensor  of  PDF.  Several 
types  of  equilibria,  all  obeying  this  symmetry,  are  found  using  numerical  continuation  software, 
including  isotropic,  prolate  and  oblate  uniaxial  and  biaxial  states.  For  equilibria  obeying  this 
symmetry,  the  PDF  reduces,  in  a  coordinate  system  specified  by  the  orthogonal  frame  of  the 
second  moment  tensor,  to  a  Boltzmann  distribution  parametrized  by  three  order  parameters  and 
material  parameters.  This  is  the  key  analytical  foundation  for  all  previous  rigorous  analyses, 
which  then  proceed  from  the  detailed  properties  of  the  Boltzmann  parametrization;  our  analysis 
likewise  begins  from  this  representation  of  the  PDF  of  the  Smoluchowski  equation  for  the 
coupled  potential.  From  a  physical  perspective  of  observable  PDFs  for  dipolar  rigid  rod 
ensembles  in  a  viscous  solvent,  the  primary  result  established  in  this  paper  is  a  proof  of 
the  numerical  observation  in  [14]:  all  stable  equilibria  for  coupled  dipolar  and  Maier-Saupe 
potentials  are  either  isotropic  or  prolate  uniaxial.  Thus,  all  oblate  uniaxial  and  biaxial  equilibria, 
which  indeed  exist,  are  necessarily  unstable.  We  remark  that  this  result  significantly  reduces 
any  numerical  computation  of  stable  equilibria,  since  one  may  posit  an  appropriate  Boltzmann 
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distribution  with  symmetry  reduced  degrees  of  freedom!  In  order  to  prove  the  main  result 
in  the  paper,  a  series  of  lemmas  is  developed.  Each  lemma  addresses  some  aspects  of  the 
probability  density  function,  which  are  essential  in  the  proof  of  the  main  theorem. 


2.  Equilibria  of  rigid,  dipolar  rod  ensembles  or  ‘extended  nematics’ 


We  denote  the  inherent  dipole  of  rigid  rods  by  a  unit  vector  m,  and  let  p{m)  denote  the 
probability  density  function  of  the  dipolar  direction  for  the  ‘extended  nematic’  defined  as 
ensembles  of  rigid  rods  with  inherent  dipoles.  The  evolution  of  the  probability  density  function 
for  the  rod  ensemble  is  governed  by  the  Smoluchowski  equation  [8]: 

dp  d  {  I  dU  dp\ 

—  =  Z) - - p  +  —\,  (1) 

dt  dm  ykgT  dm  dm J 

where  d/dm  is  the  orientational  gradient  operator  [3],  D  the  rotational  diffusivity,  kg  the 
Boltzmann  constant  and  T  the  absolute  temperature.  For  extended  nematics,  the  total  mean- 
field  interaction  potential  U (m)  consists  of  two  parts:  the  dipole-dipole  interaction  and  the 
Maier-Saupe  interaction, 

U  (m)  =  —a{m)  ■  m  —  b{m  ®m)  :m®m,  (2) 


where  a  denotes  the  strength  of  the  dipole-dipole  interaction  and  b  represents  the  strength  of 
the  Maier-Saupe  interaction.  It  should  be  pointed  out  that  the  mean-field  interaction  potential 
U (m)  is  the  part  of  the  interaction  potential  per  molecule  that  is  proportional  to  the  polymer 
concentration,  which  is  the  second  virial  coefficient  of  the  chemical  potential.  In  (2),  for 
simplicity,  we  have  normalized  the  potential  hy  kgT .  In  this  paper,  we  consider  the  case  of 
a  >  0,  which  means  energetically  a  dipole  tends  to  align  with  surrounding  dipoles. 

Equilibrium  solutions  of  the  Smoluchowski  equation  are  given  by  the  Boltzmann 
distribution  [8] 


1 

Peq(»l)  =  -exp[-f/(m)]. 


z  = 


L 


exp[— 17  {m)]dm. 


(3) 


where  Z  is  the  total  partition  function  and  S  —  {m  |  ||»i  ||  =  1 }  is  the  unit  sphere.  Notice  that  the 
Smoluchowski  equation  (1)  is  nonlinear:  the  mean- field  interaction  potential  depends  on  {m) 
and  (ffi  0  ffi)  which,  in  turn,  depends  on  the  probability  density.  As  a  result,  the  Boltzmann 
distribution  provides  an  implicit  expression  for  equilibria.  The  nonlinear  integral  equation  for 
{m)  is 


{m)  —  j  mpeqim)dm,  (4) 

while  the  nonlinear  integral  equation  for  {m  0  m)  is 

(»i0ffi)  =  J  m  <Si  mpeq{m)dm.  (5) 

Since  (m  0  m)  is  symmetric,  it  can  be  diagonalized  by  an  orthogonal  transformation. 
We  choose  a  Cartesian  coordinate  system  such  that  {m  0  m)  is  diagonal.  Furthermore,  we 
select  the  positive  direction  of  the  z-axis  such  that  if  {m)  is  not  zero,  then  {m)  has  a  positive 
component  along  the  z-axis.  This  choice  of  coordinates  removes  orientational  degeneracy  of 
all  equilibria.  Any  orthogonal  transformation  of  coordinates  yields  another  equilibrium,  so 
that  isotropic  equilibria  are  unchanged  whereas  any  anisotropic  equilibrium  corresponds  to  an 
orthogonal  group,  S^,  of  equilibria. 
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For  the  convenience  of  discussion  below,  we  introduce  notations  for  the  components  of 
{m)  and  {m  ®  m)  in  the  selected  Cartesian  coordinate  system: 

m  —  (mi,  m2,  m^). 


qi  =  (mi),  q3  >  0  if  (quqi,  qj)  +  0, 

Si  =  (m]), 

U  (m)  —  —a(q\m\  +  ^2^2  +  q^m^)  —  b(s\m\  +  S2m\  +  s^m^).  (6) 


In  this  Cartesian  coordinate  system,  an  equilibrium  is  completely  specihed  by  [qi,  Si,  i  — 
1,  2,  3}.  The  nonlinear  equations  for  [qj,  i,  }  are 


qi  —  j  m;Peq(»l)d/M, 

SijSi  —  j  mtm j Peq(m)dm. 


(7) 


(8) 


In  the  following,  we  present  all  theorems  with  respect  to  the  coordinates  above  unless  stated 
otherwise. 

First,  we  recall  the  theorem  proved  in  [14]. 


Theorem  1.  If  an  equilibrium  solution  satisfies  q^  >  0,  then  both  qi  and  q2  must  vanish. 


This  theorem  states  that  the  polarity  vector,  if  non-zero,  must  be  in  the  direction  of  one  of 
the  principal  axes  of  the  second  order  moment  tensor.  Next,  we  give  an  alternative  proof  of 
the  theorem,  which  sets  the  stage  for  subsequent  characterization  of  our  principal  focus:  the 
symmetric  subgroup  of  stable  equilibria. 


Proof  of  theorem  1.  We  prove  qi  =  0  by  contradiction.  Suppose  there  is  a  solution  satisfying 
^3  >  0  and  q\  0.  We  are  going  to  show  that  ^3  >  0  and  3:^:  0  leads  to  1/713)  >  0, 
which  contradicts  the  selection  of  the  Cartesian  coordinate  system: 

1  1  f 

—  {/M1/M3)  =  -  /  mim2exp[—U(mi,  m2,  m2)]dm 

q\  qiZ  Js 

1  f 

=  -  /  mi/M3{exp[— [/(/Ml, /M2,  m3)]  —  exp[— f/(mi,  7^2,  — m3)]}dm.  (9) 

qiZ  Jm3>0 

In  the  above  equation,  the  second  factor  of  the  integrand  is 

yi(mi,  m2,  m2)  =  exp[— {/ (mi,  m2,  m3)]  —  exp[— [/ (mi,  m2,  —m3)] 


=  2  sinh(o'g3m3)  exp(o'^imi)  exp[Q!^2m2  —  t/Ms(m)],  (10) 


where  UMs(m)  —  —b(sim\  +  S2m\  -H  s^m^)  is  the  Maier-Saupe  potential,  which  is  an  even 
function  of  mi,  m2  and  m3.  Substituting  into  (9),  we  get 


1 

—  (mimf) 

q\ 


1  f 

-  /  mim^Yiimi,  m2,  m2)dm 

q\Z  Jm': 


'  m3>0 


—  ( 

q\Z  Jmi>0,m3>0 


-  I  mim3{yi(mi,  m2,  m3)  —  yi(— mi,  m2,  m3)}d/M 


4 

Z 


/ 

t/m 


mi>0,m3>0 


mim3  sinh(Q'^3m3) 


sinh(a^imi) 

qi 


X  exp[Q!^2m2  -  UMs(m)]dm. 


(11) 


Characterization  of  stable  kinetic  equilibria  of  rigid,  dipolar  rod  ensembles 


281 


The  integrand  in  (11)  satisfies 

sinh(ar/imi) 

m\m^  sinh(a^3m3) - exp[a^2^2  ~  1^ms(»*)]  >  0 

for  mi  >0,  m3  >  0,  c/3  >  0  and  qi  ^  0.  (12) 

It  leads  to  ^(mim^)  >  0,  which  contradicts  the  selection  of  the  Cartesian  coordinate  system. 
Therefore,  we  must  have  =  0;  c/2  =  0  can  be  proved  in  a  similar  way. 

To  facilitate  the  analysis  below,  we  introduce  a  lemma. 

Lemma  1.  Let  u  —  (ui,  U2,  ■ .  ■ ,  Un)  denote  a  point  in  the  n-dimensional  space.  Consider  a 
region  A,  which  can  be  the  whole  n-dimensional  space  or  an  n-dimensional  sub-region  of  the 
n-dimensional  space  or  a  lower-dimensional  sub-region  of  the  n-dimensional  space  (such  as 
the  unit  sphere  in  the  three-dimensional  space).  Consider  a  family  of  v -dependent  probability 
density  functions  over  the  region  A: 

1 

p{u,  v)=  -  exp[u/i(M)  +  /2(m)], 

Z  =  f  exp[u/i(H)  + /2(H)]dH.  (13) 

J  A 

Let  (■)  denote  the  mean  and  var{  }  denote  the  variance,  taken  with  respect  to  the  probability 
density  p(u,  v).  Let  us  treat  v  as  an  independent  variable  and  consider  two  functions: 

gi(v)^  ifi(u)},  (14) 

g2(v)^  ifi^u)}.  (15) 

We  have  the  following 

•  lfvar{/i(M)}  >  0  for  all  V,  then  gi(ri)  is  strictly  increasing  for  all  v. 

.  lfvar{/2(H)}  >  0  for  all  V  and  g^iO)  =  0,  then  g2(v)  is  strictly  increasing  for  all  u  >  0 
and  g2(v)  is  strictly  decreasing  for  all  v  <  0. 

Remark  1.  The  two  conditions  ‘var{/i(H)}  >  0  for  all  v’  and  ‘var{/j^(M)}  >  0  for  all  v’  are 
usually  easy  to  verify.  The  only  condition  left  to  verify  is  g^iO)  =  0. 

Proof  of  lemma  1.  Differentiating  p(u,  v)  with  respect  to  v  gives 

^p(u,  v)  =  [fi(u)  -  (fi(u)}]p(u,  v).  (16) 

di) 

From  the  derivative  of  the  probability  density,  it  follows  immediately  that 

g'i(v)  =  (fi(u)[fi(u)  -  {/i(«))])  =  var{/i(«)}.  (17) 

Thus,  var{/i(H)}  >  0  for  all  v  implies  that  gi(ii)  is  strictly  increasing  for  all  v. 

To  prove  the  property  of  g2(v),  we  point  out  a  useful  fact  from  calculus. 

•  Suppose  function  f(v)  satisfies  /'(uq)  =  0  and  further  obeys  the  property  that  wherever 
f'(v*)  is  zero  f"(v*)  must  be  positive.  Then  we  have  f'{v)  >  0  for  all  v  >  Vq  and 
f'(v)  <  0  for  all  V  <  vq. 
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It  is  important  to  point  out  that  this  simple  fact  of  calculus  allows  us  to  study  the  global 
behaviour  of  a  function  from  its  local  properties;  this  observation  is  key  to  our  earlier  analysis 
of  equilibria  of  Maier-Saupe  potentials  [17].  Differentiating  g2{v)  yields 

g'(u)  =  {/2(m)[/i(«)-  {/i («))]),  (18) 

=  {/f(«)[/l(«)  -  (/l(«)>])  -  (/l"(«)[/l(M)  -  {/l(M)>]>  ■  (/l(«)) 

-{/"(«)>  ■{/l(H)[/l  («)-{/!  («)>]> 

=  {/"(«)[/"(«)  -  {/"(«)>]>  -  2{/2(h)[/i(h)  -  {/i(H)>]>{/i(M)) 

=  var{/2(H)}-2g'(u){/i(«)>.  (19) 

It  is  clear  that  function  g2{v)  satisfies  the  property  that  g^iv*)  =  0  implies  g'-^iv*)  >  0.  The 
condition  of  the  lemma  specifies  that  g^CO)  =  0-  From  the  simple  fact  of  calculus  we  pointed 
out  above  that  it  follows  that  >  0  for  all  u  >  0  and  g2(i’)  <  0  for  all  u  <  0. 

Theorem  2.  If  an  equilibrium  solution  satisfies  >  0,  then  a  and  s^  must  satisfy  as^  >  1. 

Proof  of  theorem  2.  We  prove  by  contradiction.  Suppose  asj,  ^  1.  From  theorem  1,  we 
see  that  ^3  >  0  implies  qi  —  q2  —  0.  The  parameter  ^3  has  a  fixed  value  in  the  equilibrium 
solution  pe^{m).  Here  we  rename  it  v  and  treat  it  as  an  independent  variable.  We  consider  the 
probability  density 

p{m,  v)  —  —  exp[ai>OT3  +  b{siml  +  S2m\  +  S2m\)], 

Z  =  y  exp[ai)m3  +  h(simj  +  3'2m2  +  i'3m3)]dffi. 

Note  that  p(m,  u)|„=^3  =  Psq(m).  Consider  the  function 

F{v)  =  V  —  {mf),  (20) 

where  the  average  (■)  is  taken  with  respect  to  the  probability  density  p{m,  v)  and  is  a  function 
of  V.  F{v)  satisfies  F(0)  =  0  and  Fiq^)  =  0.  We  are  going  to  show  that  as^  ^  1  implies  that 
F'(v)  >  0  for  u  e  (0,  ^3),  which  contradicts  f(0)  =  Fiq^)  —  0. 

Differentiating  with  respect  to  v  yields 

9 

— pirn,  v)  —  a(m2,  —  {m2))p{m,  v), 
dv 

^{m3)  =  a{OT3  (m3  -  {m3)))  =  a  {{mj)  -  {m3)^)  , 

F' (v)  —  I  —  a{m\)  +  a{m3)^ .  (21) 

To  show  F'(v)  >  0  for  i)  G  (0,  ^3),  we  only  need  to  show  oi{m\)  <  1  for  v  G  (0,  ^3).  Because 
of  the  assumption  o'3'3  ^  1,  we  have  Q'(m3)„=^3  ^  1.  Thus,  to  show  a  (m^)  <  Iforu  G  (0,^3), 
we  only  need  to  show  {m\)  is  strictly  increasing  for  v  >  0. 

Let  us  apply  lemma  1  with  A  —  the  unit  sphere,  fi(m)  =  am3  and  f2(in)  — 
b(s\m\  +  S2m\  +  S3m\).  When  u  =  0,  the  probability  density  p(m,  0)  is  an  even  function 
of  m3.  So  the  condition  of  part  2  of  lemma  1  is  satisfied: 

^{/i"(/M))U=o  =  {fl(m)[f3(m)  -  (/i(/M))]>|„=o 

=  oi^{m\)\y=Q  -  {ml)\„=Q{m3)\y=Q  =  0. 


(22) 
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It  follows  from  lemma  1  that  (m^)  =  ^{f\{rn))  is  strictly  increasing  for  n  >  0.  This  leads 
immediately  to  F'{v)  >  0  for  i)  e  (0,  q-^),  which  contradicts  T’(O)  =  F{q^)  —  0.  Therefore, 
a^3  >  1. 


3.  Free  energy  and  stability 


Consider  an  arbitrary  probability  density  p(m),  not  necessarily  an  equilibrium  probability 
density.  The  free  energy  of  the  probability  density  pirn)  is  given  by  [11] 


G[p]  = 


I 


p(m)  In  p(m)dm 


■  m  p{m')p(m)Am!  Am 


b 

2 


m'  ®m'  ■.  m 


0/M  p{m')p(m)Am'  Am 


f  a 

—  I  p  (m)  In  p(m)Am - (m) 

Js  2 

=  Gen,[p]  +  Gi[p]  +  G2[p], 


(Ml)  -  -(MI  I 


(23) 


where  Gent[p]  is  the  entropic  part  of  the  free  energy  and  Gi[p]  and  G2[p]  are  the  free  energy 
parts  corresponding  to  the  two  mutual  interactions:  dipole-dipole  interaction  and  the  Maier- 
Saupe  interaction,  respectively. 

In  the  above  (■)  denotes  the  mean  taken  with  respect  to  the  probability  density  p(m).  For 
the  clarity  of  analysis  below,  we  introduce  different  notations  for  means  taken  with  respect  to 
different  probability  densities. 


•  Let  (-[eq  denote  the  mean  taken  with  respect  to  the  equilibrium  probability  density. 

•  Let  (■)  denote  the  mean  taken  with  respect  to  a  general  probability  density  (for  example, 
a  perturbed  probability  density  near  the  equilibrium  probability  density). 


Now  we  prove  two  theorems  that  relate  the  second  moment,  the  dipole-dipole  interaction 
strength  (a)  and  the  stability  of  equilibria. 


Theorem  3.  If  an  equilibrium  solution  satisfies  qi  —  0  and  asi  >  1  for  i  —  I  or  i  —  2  or 
i  —  3,  then  it  is  unstable. 


Proof  of  theorem  3.  With  the  choice  of  coordinates,  only  q^  may  be  non-zero.  Without  loss 
of  generality,  we  prove  the  theorem  for  i  —  2.  To  prove  instability,  we  only  need  to  show  that 
the  free  energy  G[p]  does  not  attain  a  local  minimum  at  Pf^(m).  More  precisely,  we  show 
that  if  as2  >  1,  then  there  exists  a  perturbed  probability  density  p{m)  arbitrarily  close  to  the 
equilibrium  probability  density  Peq  (m)  such  that 

G[p(Mi)]  <  G[Peq(Ml)].  (24) 

We  consider  the  perturbed  probability  density 

P(mi)  =  (1  h-  sm2)Peq(m).  (25) 

Since  {m2>eq  =  ^2  =  0,  we  have  p(m)Am  =  1,  which  indicates  pirn)  is  a  probability 
density  function.  We  calculate  the  three  parts  of  the  free  energy  with  the  perturbed  probability 
density  p(m). 

Using  the  Taylor  expansion 

1  9 

(a  +  Ax)  ln(a  +  A2:)  —  alna  +  [lna  +  l]Ax  H - (Ax) 

2a 


(26) 
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we  have 
Gent[p]  = 


j  p  (m)  In  p(m) dm 

I  Peq(»*)lnpeq(»l)d/W  +  e  /  OT2Peq(»l)[lnpeq(»l) 

Js  Js 


+  \]dm 


+  e 


'  f  mlp^  im)  ^  ( 

Js  ^  2peq(lM) 


-dm  +  ■  ■  ■ . 


—  G‘ent[Peq]  ^  2  ^^2)eq  +  '  ‘  ■  (27) 

Because  =  0  and  q2  —  0,  the  equilibrium  probability  density  Peq(»i)  is  an  even  function  of 
m\  and  m2.  The  first  moment  of  the  perturbed  probability  density  and  Gi[p]  are 

{m)  —  {/M>eq  +  s{mm2)eq  —  (0,  0,  ^3)  +  S  (0,  (m2)eq,  0),  (28) 

Gdp]  =  -^{m)-{m) 


=  Gi[Peq]-£^^{m^)2q. 

The  second  moment  of  the  perturbed  probability  density  is 
{m®m)  —  {m®  »i)eq  +  s{mmm2)eq 


/si  0  0\  /O  0  0  \ 

=  |0  S2  o|+e(o  0  {n^3»*2)eq  1  ■ 

Vo  0  S3J  Vo  0  / 

Substituting  the  above  into  G2[p]  yields 
b 

G2[p]  =  -{m  0  m)  :  (m0m) 

=  G2[peq]  -  e^b{m3ml)l^. 

Combining  (27),  (29)  and  (31)  we  obtain 

G[p]  -  G[peq]  =  e^^(m^)eq  -  -  S^b{m3ml)l^  +  ■■■ 


(29) 


(30) 


(31) 


< -e^^£2(a£2  -  1)  +  ■  ■  ■  ■  (32) 

IfQ!£2  >  1,  then  fore  sufficiently  small  we  have  G[p]  <  G[peq],  which  implies  the  equilibrium 
solution  Peq  is  unstable. 


Theorem  4.  If  an  equilibrium  solution  satisfies  q^  >  0  and  S3  ^  S2.  then  it  is  unstable. 


Proof  of  theorem  4.  Theorem  1  shows  that  q^  >  0  implies  q^  —  q2  —  0.  It  follows  from 
theorem  2  that  q^  >  0  implies  o!J3  >  1.  Combining  these  results  with  the  condition  S2^  S3, 
we  obtain  q2  —  0  and  as2  0CS3  >  1.  Theorem  3  implies  that  the  equilibrium  solution  is 
unstable. 
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4.  All  stable  solutions  are  axisymmetric 

In  this  section,  we  show  that  all  stable  equilibrium  solutions  must  be  prolate  uniaxial 
(axisymmetric).  Note  that  the  stability  requirement  is  essential  in  establishing  the  axisymmetry. 
The  existence  of  unstable  non-axisymmetric  equilibrium  solutions  is  presented  in  the  appendix. 

We  first  recall  the  result  proved  for  pure  (non-dipolar)  nematic  rod  ensembles  where 
{m)  =  0  [11,19]. 

Theorems.  For  pure  nematics,  the  stable  equilibria  satisfy  either  s\  —  S2  —  so,  (isotropic 
phase)  or  ^3  >  sj  =  S2  (prolate  phase)  in  the  selected  coordinate  system. 

Below  for  the  case  of  extended  nematic  equilibria  in  which  {m)  may  be  non-zero,  we 
are  going  to  show  that  if  {m)  ^  0  (i.e.  qi  —  q2  —  Q  and  ^3  >  0  by  the  selection  of  the 
coordinate  system  and  by  the  result  of  theorem  1),  then  a  stable  equilibrium  solution  must  be 
uniaxial.  Furthermore,  the  axis  of  symmetry  must  be  the  major  director  (i.e.  the  eigenvector 
of  the  second  moment  corresponding  to  the  largest  eigenvalue).  That  is,  {m\)  —  (m\)  <  {m\). 

Before  we  go  into  the  details  of  proving  the  axisymmetry,  let  us  describe  the  general 
approach  we  are  going  to  use  in  the  proof.  A  key  step  in  the  proof  for  the  case  of  pure  (non- 
dipolar)  nematics  is  that  we  renamed  as  n  and  treated  v  as  an  independent  variable  [19].  For 
pure  nematics,  we  parametrized  the  probability  function  by  introducing 

n  =  Si  -  (mf)y=o  (33) 

and  proved  that  hj  =  —  {mj)  has  the  property 

{hfyhi  -  h2))  <  0  ifri>r2.  (34) 

Then  the  function 

F(v)  =  r2({m\)  -  {mi)„=o)  -  ri({ml)  -  {m\)y=o) 

satisfies  F(0)  —  F(b)  —  0,  which  will  lead  to  a  contradiction  if  the  axisymmetry  is  violated. 
For  the  extended  nematics,  the  definition  of  r,  is  non-trivial.  Note  that  for  a  definition  of  r, 
to  work  in  the  proof,  it  has  to  preserve  property  (34).  For  the  extended  nematics,  definition 
(33)  does  not  preserve  property  (34).  Lemmas  2-5  presented  below  are  the  intermediate  steps 
leading  to  a  successful  definition  of  r,-  for  the  extended  nematics.  Then  finally  in  theorem  6, 
the  results  of  these  lemmas  are  used  to  prove  the  axisymmetry. 

Lemma  2.  Consider  a  family  of  r -dependent  probability  density  functions  of  9 
1  2 

p(0,  r)  =  —  exp[flir  cos20  cos0], 

p27l 

Z=  /  exp[flir^  cos20  H- cos0]d0.  (35) 

Jo 

Let 


g{r,ai,a4)  =  (cos  26),  (36) 

where  the  average  is  taken  with  respect  to  the  probability  density  (35).  If  ai  >  0,  then 
g{r,  ai ,  04)  is  a  strictly  increasing  function  of  r  for  r  >  0. 

Note  two  aspects  of  this  lemma.  First,  there  is  no  condition  imposed  on  04.  Second,  a 
special  case  of  lemma  2  with  04  =  0  was  proved  in  [5]  and  was  used  in  [19]  to  prove  that  all 
equilibria  of  pure  (non-dipolar)  nematic  polymers  are  axisymmetric. 


286 


H  Zhou  et  al 


Proof  of  lemma  2.  To  prove  lemma  2,  let  us  replace  by  Vi,  replace  r  by  V2  and  treat  Ui  and 
V2  as  two  independent  variables.  ai  and  <24  are  treated  as  parameters.  The  probability  density 
becomes 

1 

p{9,  Vi,  V2)  =  —  exp[flii>i  cos 20  +04112  cos0], 

Z—  I  exp[oii;i  cos20 +a4V2Cos0]d0.  (37) 

Jo 

Consider  the  function 

gsivu  V2)  =  (cos 20).  (38) 

We  first  fix  V2  and  treat  Vi  as  a  variable.  Let  us  apply  lemma  1  with  A  =  [0,  27r], 
fi{6)  =  aicos20  and  /2(0)  =  04^2 cos 0.  Clearly,  var(/i)  >  0.  It  follows  that 
{/i(^))  =  cii{cos20)  =  01^3(111,  V2)  is  strictly  increasing  for  all  vi,  which  means  g3(i>i,  V2) 
is  a  strictly  increasing  function  of  Vi  for  all  Ui. 

If  04  =  0,  then  g3(ni ,  +2)  does  not  depend  on  V2.  Suppose  04  ^  0.  We  fix  vi  and  treat  112  as 
a  variable.  Let  us  apply  lemma  1  with  A  =  [0,  27r],  fi  (0)  =  04  cos  0  and  /2(0)  =  oi  Vi  cos  20. 
When  V2  =  0,  expressing  cos  26  —  2  cos^  0  —  1,  we  see  that  the  probability  density  at  i;2  =  0 
is  an  even  function  of  cos0  which  implies  (cos^  6)\y.^=o  —  0  for  any  odd  integer  k.  So  the 
condition  of  part  2  of  lemma  1  is  satished:  var(/j^)  >  0  and 

^(fHo))\v,=o  =  (fhmm  -  {/i(0)>])u.=o 
av2 

=  fl4{cos^0)|„2=o  -  a4{cos^0)|„2=o{cos0)|„j=o  =  0.  (39) 

It  follows  that  if ^(6))  =  alicos^  0)  =  Q2g3(i'i+2)+i  strictly  increasing  for  all  V2  >  0,  which 
means  §3(^1,  V2)  is  a  strictly  increasing  function  of  V2  for  all  V2  >  0  if  04  ^  0.  If  04  =  0, 
g3(fi,  V2)  does  not  depend  on  +2- 

By  definition,  function  ^3  is  related  to  g  by 

g{r,aua4)  =  g3(r^,r).  (40) 

Since  ^3(1^1,  V2)  is  a  strictly  increasing  function  of  Vi  for  all  Vi  and  is  either  independent  of 
V2  (if  04  =  0)  or  a  strictly  increasing  function  of  V2  for  1)2  >  0  (if  04  ^  0),  we  conclude  that 
g(r,  a\ ,  04)  is  a  strictly  increasing  function  of  r  for  r  >  0. 


Lemma  3.  Consider  a  family  of  v -dependent  probability  density  functions  of  m: 
p(m,  n)  =  —  exp[c4OT3  +  Cim^  +  v(c2m\  +  C3m3)], 

Z  =  y  exp[c4m3  +  cim3  +  i;(c2nr2  +  C3m3)]d»i.  (41) 

Let  (■)  denote  the  average  taken  with  respect  to  p(m,  11).  Let  h /  =  mj  —  (mj).  hj  has  the 
properties 

h  \  +  /?2  +  ^3  —  0, 

{/zi>  =  (h2)  =  (h)  =  0.  (42) 

Suppose  (ci,  C2,  C3,  C4)  satisfy  the  condition  ci  +  VC3  >  VC2  >  0  for  all  v  >  0.  For  v  >  0, 
we  have 


(hiiM  -  h2))  <  0, 
{h2ih3  -  hi))  <  0, 
{h3ih2  -  hi))  <  0. 


(43) 
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Remark.  (hi{h^  —  h2))  is  the  correlation  of  hi  and  —  h2). 


Proof  of  lemma  3.  Below  we  present  the  proof  of  {hiih^  —  h2))  <  0.  {/!2(/*3  —  hi))  <  0  and 
(/!3(/z2  —  hi))  <0  can  he  proved  in  a  similar  way. 

Expressing  in  terms  of  nii,  m2,  we  have 

{hi(h2  -  h2))  =  {{m\  -  {m\)){ml  -  mj)).  (44) 

We  select  the  x-axis  (the  axis  associated  with  mi)  as  the  direction  of  the  pole  of  a  spherical 
coordinate  system  (0,  9).  In  the  spherical  coordinate  system  (0,  9),  we  have 

mi=cos0,  m2  =  sin(/>  sin0,  m3  =  sin 0 cost?, 

m\  —  cos^  (j), 

m\  —  m\  —  sin^  <p  cos  29, 

C4m2  +  (ci  +  vc2)m\  +  vc2m\  —  C4  sin  <p  cos  9  +ai  sin^  (p  cos  29  +  02  sin^  <P , 

(45) 


where  the  coefficients  ai  and  02  are 

Cl  +  r(c3  -  C2) 
«i  = - ;; - 


>  0, 


fl2  = 


Cl  +  V(C3  +  C2) 


(46) 


The  probability  density  function  (41)  becomes 

p(0,  9)  —  —  exp[fli  sin^  (p  cos  29  +  C4  sin 0  cost?]  exp[a2  sin^  (p]. 


Z  = 


ff 


exp[ai  sin^  (p  cos  29  +  C4  sin  (p  cos  0]d0  exp[fl2  sin^  (p]  sin  (p  dcp. 

(47) 


We  compare  exp[ai  sin^0cos20  +  C4  sini/;  cos  0]  with  the  probability  density  in  lemma  2 
and  identify  (C4,  sin0)  with  (04,  r)  in  lemma  2.  It  follows  that  p{(p,  9)/  p((p,  9)d9  is  a 

probability  density  function  of  9  and  has  the  same  form  as  the  one  in  lemma  2.  Applying 
lemma  2  gives 


p2n 

L  ' 


cos20p(0,  0)d0  =  g(sin0,  ai,  C4)  /  p{(p,9)d9. 


pin 

Jo 


(48) 


where  g{r,  ai ,  04)  is  defined  in  (36),  and  for  ai  >  0,  gir,  ai,  04)  is  a  strictly  increasing  function 
of  r  for  r  >  0. 

Let  cos^  (po  —  (cos^  cp).  Writing  (/zi(/z3  —  02))  in  spherical  coordinates,  we  have 
{hi(h2-h2))  =  {(m\  -  {m\)){ml  -  m\)) 


—  I  (cos"^  (p  —  cos^  (po)  sin^  (p 

Jo 


j  cos  20p(0,  0)d0^  sin  0  d(/> 


(cos^  (p  —  cos^  (po)  sin^  (pgisincp,  ai,  C4) 


pin 

Jo 


p((p,  9)d9  I  s'mcpdip 


—  ((cos^  (p  —  cos^  (po)  sin^  (pg{sm(p,  ai,  C4)) 

—  ((cos^  (p  —  cos^  (po)[sm^  (pgisiiKp,  ai,  C4)  —  sin^  (pog{&m(po,  ai,  C4)])  <  0. 
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In  the  above  we  have  used  the  fact  that  sin^  0g(sin0,  a\,  C4)  is  a  strictly  increasing  function 
of  sin  (p  and  that  cos^  0=1  —  sin^  0  is  a  strictly  decreasing  function  of  sin  0. 

Lemma  4.  Consider  a  family  of  v -dependent  probability  density  functions  of  m: 
p(m,  v)  —  —  exp[c4m3  +  ciml  +  v{c2m\  +  C2,m\)], 

Z  =  y  exp[c4m3  +  cim3  +  nfewij  +  C3m3)]d»i. 

Let  (■)  denote  the  average  taken  with  respect  to  p{m,  v).  Suppose  (ci,  C2,  C3,  C4)  satisfy  the 
condition  ci  +  t;c3  >  VC2  >  0  for  all  u  >  0  and  C4  ^  0.  Then  we  have 

L  {m\)y=Q  is  a  strictly  increasing  function  of  ci  for  all  ci; 

2.  {ml)^=Q  >  5  >  {m\)v=o  =  (m^)„=o; 

3.  for  n  >  0,  {mj)  -  {ml)y=o  >  {ml)  -  {ml)y=o. 


Remark.  Note  that  condition  c\  +  vcs  >  VC2  >  0  for  all  u  >  0  Implies  that  ci  ^  0, 
C3  ^  C2  >  0. 


Proof  of  lemma  4. 

Proof  of  item  1.  To  prove  that  {ml)y=Q  is  a  strictly  increasing  function  of  ci  for  all  ci,  we 
treat  ci  as  the  independent  variable  and  hx  i;  =  0  and  other  parameters.  Applying  lemma  1 
with  A  =  unit  sphere,  f\{m)  —  ml,  and  f2im)  —  C4OT3,  we  arrive  at  that  {ml)y=o  is  a  strictly 
increasing  function  of  ci  for  all  ci. 

Proof  of  item  2.  For  ci  ^  0,  using  item  1,  we  have 

{ml)y=c^=o  >  (m3)„=ci=c4=o-  (49) 

Next  we  treat  C4  as  the  independent  variable  and  hx  v  =  0  and  other  parameters.  We 
apply  lemma  1  with  A  =  unit  sphere,  f\(m)  —  m3  and  /2(»i)  =  cim^.  It  is  straightforward 
to  verify  that  the  condition  of  lemma  1  is  satished.  Thus,  we  obtain  that  {ml)y=o  is  a  strictly 
increasing  function  of  C4  for  all  C4  >  0  and  a  strictly  decreasing  function  of  C4  for  all  C4  <  0. 
That  is,  {ml)y=o  attains  the  global  minimum  at  C4  =  0.  Thus,  for  C4  ^  0,  we  have 

{ml)y=o  >  {ml)y=c^=Q.  (50) 

When  Cl  =  C4  =  n  =  0,  the  probability  density  gives  a  uniform  distribution,  which  leads  to 

{m  f)  y—c^—C^  —  Q  {mf)y—y^—y^  —  Q  {m  f)  y  —  y^—y^  —  Q  3.  (51) 

Combining  (49),  (50)  and  (51),  we  obtain 

{mf)y=(2  >  {mf)y  =  y^  =  Q  ^  {m  f)  y  =  y^=y,^  =  Q  —  3.  (52) 

Using  the  fact  m\+ml+ml  =  1  and  the  fact  that  at  n  =  0  the  probability  density  is  independent 
of  mi  and  m2,  we  have 

{m\)y=Q  =  (m2)„=o  =  5(1  -  {ml)y=o)  <  5-  (53) 

Proof  of  item  3.  To  prove  item  3,  we  only  need  to  show  that  for  u  >  0 
{m3)  -  {ml)  >  {ml)y=o  -  {ml)y=o- 


(54) 
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That  is,  (m3)  —  (m^)  is  a  strictly  increasing  function  of  v.  Recall  that  hj  is  defined  as 
hj  —  m^j  —  {m^j)  and  satisfies  (hj)  —  0.  Taking  the  derivative  with  respect  to  v,  we  obtain 

^p{m,  v)  =  [C3(m3  -  (m3))  +  C2(m2  -  (m2})]p(m,  v) 
av 

=  [cshs  +C2h2\p(m,  n), 

^  (m?  -  m\)  =  {{ml  -  m\){c2h2  +  Ci/ii)) 
dn 

=  {{h  -  /l2)(C3/l3  +  C2h2)) 

=  {(h3  -  h2)  +  h2)  +  -  112)^1  . 

Using  /z3  +  /z2  =  —hi  and  using  the  result  of  lemma  3,  we  get 

d  0  0  Co  C'X  —  Co  0 

{ml  -  ml)  =  -±^{h,{h2  -  h2))  +  ^— ^((/13  -  h2)^)  >  0.  (55) 

dn  2  2 

Here  we  have  used  (c3  +  C2)/2  >  0,  (c3  —  C2)/2  ^  0,  {hiih^  —  h2))  <  0  and  {{h^  —  h2)^)  >  0. 
Lemma  5.  Consider  the  equilibrium  probability  density: 

Peq(/«)  =  ^  exp[Q!53m3  +  b{{s2  -  Si)ml  +  (53  -  ^Om^)], 

Z  =  y  exp[ci'^3m3  +  ii((i2  —  ■?i)m2  +  (^3  —  ii)m3)]d/M.  (56) 

We  first  define  functions  ci  (ct),  C2((t),  03(0-)  and  04(0)  as  follows: 

3cr  -  1 

ci(cr)  =  b — - — , 

C2(cr)  =  ^2  -  ^1, 

3ct  -  1 

C3(cr)  =  ^3  -  - - — , 

C4{a)  =  aq2,  (57) 

where  q^  and  Sj  {i  =  1,  2,  3)  are  the  first  and  second  moments  of  the  equilibrium  solution, 
respectively.  We  construct  a  family  of  (a,  i;)-dependent  probability  density  functions  of  the 
form  (41): 

p{m,  cj,v)  —  —  exp[c4(CT)m3  +  ci(cr)m3  +  n(c2(cr)m2  +  C3(cr)m3)], 

Z  =  y  exp[c4(CT)m3  +  ci(cr)m3  +  i;(c2(cr)m2  +  C3(CT)m3)]d/M.  (58) 

It  is  straightforward  to  verify  that 

p  (m ,  (T,  n)  I  „=z,  =  Peq  (m) .  (59) 

Let  (■)„  be  the  mean  taken  with  respect  to  probability  density  p{m,cT,  v).  We  define  a  mapping 
p  =  H(a)  as 

q  ^  H{a)  =  {ml)„\y=o-  (60) 

Using  the  value  of  q,  we  set  parameters  ci  (q),  C2{q),  C2,{q)  and  C4(?7)  and  construct  probability 
density  p{m,  q,  n)  by  replacing  Cj{a)  in  (58)  with  cjiq). 
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Suppose  ^3  >  0.  If  CT  satisfies  the  conditions 

1.  cr  ^  i  and  H{<j)  >  a  (i.e.  rj  >  a), 

2.  ci(ct)  +  ucsCct)  >  vc2(cr)  >  0  for  all  i;  >  0, 

then  rj  =  H (a)  has  the  properties 

1.  p  >  i  and  H(ri)  >  rj, 

2.  ci(rj)  +  vc-iirj)  >  vc2(rj)  >  0  for  all  u  >  0. 

Remark.  Mapping  rj  =  H(a)  preserves  the  property  ci  +  vcj,  >  VC2  >  0  for  all  i>  >  0.  As 
long  as  this  condition  is  satisfied,  lemmas  3  and  4  can  he  applied.  As  we  will  see,  lemmas  3 
and  4  will  play  an  essential  role  in  the  proof  of  the  main  result  that  all  stable  equilibrium 
solutions  are  either  isotropic  or  prolate  uniaxial. 

Proof  of  lemma  5. 

Proof  of  item  1.  From  item  1  of  the  condition:  cr  ^  |  and  rj  —  H(a)  >  cr,  we  have  p  >  |  and 
3n  -  1  3ct  -  1 

ci(rj)^b — - —  >b — - — =ci(ct)^0.  (61) 

Note  that  C4(ct)  =  aq^  is  independent  of  a  and  that  {m^)y=o  is  independent  of  C2  and  C3.  In 
other  words,  depends  on  ct  only  through  its  dependence  on  ci.  Result  1  of  lemma  4 

shows  that  {m^)„=Q  is  a  strictly  increasing  function  of  ci.  Thus,  ci{rj)  >  ci(cr)  gives  us 

H(rj)  =  {m3)J„=o  >  {ml)a\v=o  =  H((j)  =  rj.  (62) 

Proof  of  item  2.  Now  we  prove  Ci(p)  +  vc^frj)  >  vc2(rj)  >  0  for  u  >  0.  Notice  that 
C2(cr)  =  S2  —  ^1  is  independent  of  ct.  So  vc2(rj)  >  0  for  u  >  0  follows  directly  from 
uc2(cr)  >  0  for  i;  >  0.  We  have  already  obtained  ci{rj)  >  0.  Hence,  we  only  need  to  prove 
C3('?)  -  C2(rj)  >  0. 

Item  2  of  the  condition,  ci(ct)  +  vc^ia)  >  vc2(cr)  >  0  for  all  u  >  0,  allows  us  to  apply 
lemma  4  with  probability  density  p{m,  a,  v).  Using  result  3  of  lemma  4  at  it  =  we  have 

1 1;=0  — 

2  2  1  —  p 

2  2  3p  —  1 

^3  -  {mf)a\v=Q  -  {S2  -  (»J2)aU=o)  =  ^3-^2 - ^ -  >  0, 

which  immediately  leads  to 

-  C2(v)  =  i'S  -  ‘52  -  >  0-  (63) 

Theorem  6.  If  a  stable  equilibrium  solution  satisfies  q^  >  0,  then  s^  >  si  —  S2.  Thus,  all 
stable  anisotropic  equilibria  obey  prolate,  uniaxial  symmetry. 

Proof  of  theorem  6.  We  prove  the  conclusion  in  several  steps. 

Theorem  4  indicates  that  if  an  equilibrium  solution  with  ^3  >  0  is  stable,  then  we  have 
^3  >  i’l  and  53  >  S2.  We  need  to  show  ii  =  S2.  We  prove  it  by  contradiction.  Suppose  53  >  ■Ji- 

Step  1.  Let  ctq  =  |.  With  the  assumption  S2  >  5 1,  it  is  straightforward  to  verify  that 
ci(cro)  +  uc3(o-o)  >  vc2(o'o)  >  0  for  u  >  0.  This  allows  us  to  apply  lemma  4.  Result  2  of 
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lemma  4  implies  that  H(ao)  —  (m^)crfj\v=o  >  \  —  oq-  Thus,  the  two  conditions  of  lemma  5 
are  satisfied  for  (Tq  =  f  Starting  from  oq  —  I,  we  construct  a  sequence  {cr„}; 


Gk  =  H{ak-i). 


(64) 


Lemma  5  guarantees  that  for  each  new  the  two  conditions  of  lemma  5  are  satisfied  and  that 
Oji  >  Ok-i-  Since  is  the  average  of  m\,  we  have  ^  1.  So  the  sequence  {aic}  is  strictly 
increasing  and  bounded.  Hence  the  sequence  {ak}  has  a  well-defined  limit  as  A:  ^  oo.  Let 
a  —  limi:^oo  <^k-  Lemma  5  implies  that  a  >  |.  Let  us  introduce 


We  have 


r-i  =  s^- 

r2  =  S2 

ri=si- 

Ci(ct)  =  b- 


CT; 

1  —  O’ 

'  2  ^ 
1  —  cr 


3cr  -  1 


>  0, 


C2(o-)  =  ^2  -  .Jl  =  ^2  -  ^1  >0, 
3a  -  1 


C3(a)  =  53  -  51  - 


rs  -  ri . 


Since  the  mapping  —  H(ak-\)  preserves  the  property  ci  +  vcj  >  VC2  >  0  for  u  >  0,  in  the 
limit  of  ^  ^  oo,  we  have  C3(a)  ^  C2(a).  We  already  know  that  ci(a)  >  0  and  C2(a)  >  0. 
Combining  these  results,  we  have 

ri+r2  +  r2  =  0, 

r3  ^  r2  >  ri , 

r3  >0,  ri  <  0, 

and  ci(a),  C2(a),  C3(a)  satisfy 

ci(a)  -I-  i>C3(a)  >  i>C2(a)  >  0,  for  u  >  0. 

Thus,  lemmas  3  and  4  can  be  applied  to  the  limit  probability  density,  which  has  the  form 


1  3a  - 

p(m,  a,  v)  —  —  exp[ag3m3  +  b — — 


1 


-mj  -H  viir2  -  ri)ml  +  (r^  -  ri)ml)], 


-L 


Z  —  l  enjp[aq2,m2,  +  b- 


3a  -  1 


-m\  -H  u((r2  —  ri)m2  +  (r^,  —  ri)ml)]dm. 


(65) 


The  limit  probability  density  p(m,  a,  v)  satisfies 


{ml)\^=o  =  a, 

{^l)\v=0  —  (^2)  1 11=0  — 


1  -  {mj)\y=o 


1  —  a 


(^in2)\y=b  (m2)|i;=0  —  ^2 


{kkli)\v=b  (mi)lu— 0  —  5l 


2 

1  —  a 
2 

1  —  a 


=  r2, 


n- 


Here  for  simplicity  we  write  (-ja-  simply  as 
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Step  2.  Consider  the  function 

F(v)  =  r2{{m\)  -  (mi)„=o)  -  -  {m\)y=Q).  (66) 

F{v)  satisfies  Fijd)  —  F(b)  —  0.  Using  the  supposition  S2  >  ^i,  we  are  going  to  show  that 
F'{v)  <  0  for  n  >  0,  which  contradicts  F(0)  =  F{b)  =  0.  Let  ,■  =  m]  —  {rrP'-).  (hi,  h2,  hj) 
satisfies 

hi  +  h2  +  ^3  —  0, 

{hi)  =  {h2)  =  (fis)  =  0.  (67) 

Using  the  result  of  lemma  3,  we  have 
{hi{h2-h2))  <0, 

{h2{h2-hi))  <0, 

{h2{h2-hi))  <0, 

^p{m,  (T,  v)  =  [(r3  -  ri)iml  -  {m]))  +  (ri  -  ri)im\  -  {m\))]p(m,  v) 
av 

=  {r^h^  +  r2h2  +  rihi)p{m,  v) 

=  [('•3  -  r2)h2  +  {ri  -  r2)hi]p(m,  v),  (68) 

F'(v)  =  r2(mi[(r3  -  r2)h'}  +  (n  -  r2)fii]>  -  ri{ml[(r2  -  r2)h2  +  (n  -  r2)hi]) 

=  {(.r2hi  -  rih2)[{r2  -  r2)h2  +  (n  -  r2)hi\).  (69) 

Using  ri+r2  —  —r^,  we  have 

r2hi  -  rih2  =  \{r2{h2  -  hi)  +  (r2  -  ri){hi  +/J2)]-  (70) 

On  the  other  hand,  using  r2  +  ^3  =  —ri  and  adding  ri(hi  +  h2  +  h^)  =  0  to  it,  we  have 
r2hi  -  rih2  =  \[{r2  -  r2)hi  +  {r2  +  r2,)hi  -  2rih2\ 

=  \[{r2  -  r2,)hi  -  rihi  -  2rih2+ri{hi  +  h2  +  h^)] 

=  2[(o  -  r2)(h2  +  hi)  +  ri(h2  -  h2)]. 

Substituting  these  two  expressions  for  ^2/11  —  rih2  into  F'(v),  we  obtain 
F'iv)  =  i([r3(fi2  -  hi)  +  {r2  -  ri)(fii  + /Z2)](r3  -  r2)hi) 

+  5([(0  -  r2)ih2+hi)+ri{hi  -  h2)](ri  -  r2)hi) 

=  20(0  -  r2){hi{h2  -  hi))  +  ^{r2  -  ri)(ri  -  r2){hi(h2  +  hi)) 

-  5(^2  -  ''i)(o  -  r2){hi(h2  +  hi))  +  l(-ri)(r2  -  ri){hi{hi  -  hi)) 

=  2^3(''3  -  r2)(hiih2  -  hi))  +  \{r2  -  n)(o  -  r2){h2(hi  -  hi)) 

+  ii-ri)(r2-ri){hi{hi-h2)).  (71) 

In  the  above  expression,  all  coefficients  are  non-negative  and  all  correlations  are  negative.  The 
coefficient  of  the  third  term  (— ri)(r2  —  ^i)  is  strictly  positive.  So  the  sum  must  be  strictly 
negative.  That  is,  F'{v)  <  0  for  v  >  0,  which  contradicts  F(0)  =  F(b)  =  0.  Therefore,  we 
conclude  i'l  =  S2- 
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For  an  axisymmetric  equilibrium  with  the  z-axis  as  the  axis  of  symmetry,  we  have 

1  -  53 

So  a  stable  equilibrium  state  of  the  extended  nematic  polymers  is  completely  specified  by  the 
values  of  and  53.  The  nonlinear  integral  equation  for  and  53  is 

93 


=  j  m^peq(m)dm, 
53  =  y  mlpeq(m)dm, 


(72) 


where  the  equilibrium  probability  density  is  given  by 


1 

Peq(»l)  =  -  exp 


aq^m^  +  b 


353  -  1 


Z  = 


I 


exp 


aqo,m-i  +  b 


353  -  1 


dm. 


(73) 


Numerical  solutions  of  equation  (72)  have  been  carried  out  in  [14].  The  solution  structure, 
affected  by  two  parameters  a  and  b,  appears  to  be  significantly  more  complicated  than  that  of 
pure  (non-dipolar)  nematics.  The  theoretical  study  of  the  solution  structure  of  equation  (72) 
will  be  the  subject  of  a  future  work. 


5.  Conclusions 

The  stable  equilibria  of  rigid,  dipolar  rod  ensembles  (so-called  extended  nematics)  are  shown 
to  be  either  isotropic  or  prolate  uniaxial.  Furthermore,  the  distinguished  axis  of  symmetry  of 
stable  anisotropic  equilibria  is  identical  to  the  first  moment  of  the  PDF  (the  polarity  vector) 
and  the  major  director  of  the  second  moment  of  the  PDF  (eigenvector  associated  with  the 
largest  principal  value).  The  property  of  axisymmetry  greatly  simplifies  any  procedure  for 
obtaining  physically  observable  equilibria,  whether  it  is  numerical  or  some  form  of  asymptotic 
approximation  scheme  involving  closure  models.  We  anticipate  the  analytical  approach 
developed  in  this  paper  can  be  extended  to  Smoluchowski  equations  with  more  complex 
potentials  or  with  external  fields. 
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Appendix.  Existence  of  unstable  non-axisymmetric  equilibrium  solutions 

In  this  appendix,  we  show  that  the  extended  nematic  polymers  do  admit  non-axisymmetric 
equilibrium  states  if  the  stability  requirement  is  removed.  Specifically,  we  only  need  to  show 
the  existence  of  a  non-axisymmetric  equilibrium  state  for  one  set  of  {b,a). 
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Theorem  1  shows  that  if  ^3  ^  0  then  we  must  have  qy  —  q2  —  Q.  Using  this  result,  we 
simplify  the  nonlinear  equations  (7)  and  (8)  to 


^3  =  y  m2Peqim)dm, 

(74) 

Si  —  y  mj Peqim)dm,  i  —1,2,3 

(75) 

where  the  equilibrium  probability  density  is  given  by 

Psq{m)  =  —  exp[Q!^3m3  +b{s\m\  +  S2m\+ 


Z  =  y  &x^{aq2m2  +  b{sim\+ S2m\+ S2m^)'\dm.  (76) 

Note  that  the  equilibrium  probability  density  Peqim)  depends  on  a  and  q^  only  through  its 
dependence  on  the  combination  X  =  aq^.  Once  the  value  of  1  =  aq^  is  given,  peq(m)  no 
longer  depends  on  the  individual  value  of  a  or  q^.  The  key  in  constructing  non-axisymmetric 
equilibrium  solutions  is  to  treat  X  =  aq^asa  parameter  and  to  satisfy  equation  (74)  by  adjusting 
the  value  of  a.  Mathematically,  with  the  introduction  of  k  =  aq^  as  the  new  parameter  and 
with  a  as  the  new  unknown,  equation  (74)  becomes 


j  m2Peq(m)dm. 


(77) 


Using  the  constraint  i'l  +  53  +  ^3  =  1,  we  see  that  only  the  first  two  component  equations  in 
(75)  are  needed;  the  third  component  equation  follows  automatically  from  the  first  two.  We 
keep  i'l  and  S2  as  the  unknowns  and  write  ^3  =  1  —  —  S2-  Below  is  our  new  formulation: 


•  we  treat  5i,  3'2  and  a  as  unknowns; 

•  we  treat  b  and  X  as  parameters; 

•  we  have  three  equations:  (77)  plus  the  first  two  components  of  (75). 

The  equilibrium  probability  density  has  the  form 

Peqim',  i'l,  S2',  b,X)  —  —  exp[km3  +  b(s\m\  +  S2m\  +  (1  —  i'l  —  S2)m\)'\, 

Z  =  y  ex\i[Xm2  +  b{s\m\  +  S2m\  + —  Si  —  S2)m\)]dm.  (78) 

If  1  =  0,  then  both  sides  of  (77)  are  identically  zero  and  equation  (77)  is  automatically  satisfied. 
If  1  ^  0,  the  equilibrium  probability  density  can  be  written  as 


Peqim)  —  esrpiXmo,)  x  (an  even  function  of  m). 

Hence,  the  right  side  of  (77)  is  non-zero  and  is  of  the  same  sign  as  k.  It  follows  that  equation 
(77)  can  always  be  satisfied  by  selecting  a  suitable  positive  value  for  a.  Therefore,  in  our 
new  formulation  with  X  =  aq^  as  a  parameter  and  a  as  an  unknown,  we  gain  the  freedom  of 
being  able  to  solve  for  a  separately  from  equation  (77)  after  (^i ,  ^2)  has  been  determined.  The 
equation  for  (^1,^2)  is 

si—j  m\peqim)Am,  S2  —  \  m\peqim)Am.  (79) 

J  s  Js 

Note  that  in  the  new  formulation,  a  does  not  appear  in  the  equation  for  (ji,  53)  (79).  This 
feature  of  the  new  formulation  is  very  useful  for  our  purpose  of  demonstrating  the  existence 
of  non-axisymmetric  equilibrium  solutions.  Before  we  go  into  details,  let  us  outline  the  steps 
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we  will  take  in  constructing  a  non-axisymmetric  equilibrium  solution. 

•  We  first  solve  equation  (79)  for  k  =  0.  Specifically,  we  take  b  —  %  and  pick  the  prolate 
solution  ^2  whose  major  director  is  parallel  to  the  y-axis. 

•  For  general  k,  we  write  equation  (79)  as  a  variational  problem  by  introducing  energy 
function  Gi(yi,  S2,  k).  We  prove  that  (y{*'\  (the  solution  we  pick  for  A.  =  0)  is  an 
isolated  global  minimizer  of  Gi  (ji ,  ^2,  0). 

•  Since  (yj*'\  is  an  isolated  global  minimizer  of  Gi(si,  J2,  0),  we  can  find  a  small 
rectangle  centred  around  (s}°\  ^2°^)  such  that  the  minimum  of  Gi(ji,  S2,  0)  on  the  boundary 
of  rectangle  is  strictly  larger  than  Gi(yj°\  0). 

•  We  show  that  9Gi  (ii  ,S2,k) /dk  is  uniformly  bounded  by  1  for  all  values  of  (ii ,  S2,  k).  So 
we  can  find  a  small  positive  Ai  such  that  the  minimum  of  Gi (ii ,  52,  ■^i)  on  the  boundary 
of  rectangle  is  strictly  larger  than  Gi(s}°\  52°^  ^i),  which  implies  that  as  a  function  of 
(5i,  S2),  Gi(si,  52,  ■^1)  attains  a  minimum  inside  the  rectangle.  Therefore,  for  A.  =  Ai 
equation  (79)  has  a  solution  near  (5^°^  5^°^)  (the  solution  we  pick  for  1  =  0). 

•  Finally,  we  show  that  the  solution  we  constructed  in  this  way  for  A  =  Ai  is  indeed  non- 
axisymmetric. 


We  start  with  the  case  of  A.  =  0.  In  this  case,  equation  (79)  with  equilibrium  probability 
density  given  in  (78)  reduces  to  the  case  of  pure  (non-dipolar)  nematic  polymers.  From  the 
theoretical  results  of  pure  nematic  polymers  [11, 15, 19],  we  know  that  for  b  >  7.5  there  is 
only  one  family  of  rotation-equivalent  prolate  equilibrium  states.  In  the  analysis  below  we  fix 
b  —  We  select  the  prolate  equilibrium  state  whose  major  director  is  parallel  to  the  y-axis. 
Thus,  for  A  =  0,  a  solution  of  equation  (79)  is  given  by 


.(0) 


1  1 


(0) 


(0)  1  ^  (0)  ^ 
5)  =  -  H-  -5^w  >  _ 

233  3 


(80) 


where  5*^°'  =  —  j)  >  0  is  the  order  parameter  of  the  prolate  solution.  Here  we  use 

the  notations  (5{*'\  to  show  explicitly  that  the  solution  corresponds  to  the  case  of  A  =  0. 
In  order  to  construct  a  non-axisymmetric  solution,  we  introduce  two  energy  functions  for 
probability  density  of  the  form  (78),  which  is  completely  specified  by  (51 ,52:  A).  Recall  that  b 
has  been  fixed  at  (7  =  8.  We  first  rewrite  equation  (79)  as  a  variational  problem  by  introducing 
the  energy  function 


Gi(5i,  52;  A) 


b 

2 


(5j  -t  52  "t"  (1  —  5l  —  52)  )  —  In  Z, 


(81) 


where  the  partition  function  Z  is  given  in  (78).  Here  5i  and  52  are  treated  as  independent 
variables  while  b  and  A  are  treated  as  parameters.  It  is  straightforward  to  verify  that  a  critical 
point  of  energy  function  Gi  is  equivalent  to  a  solution  of  equation  (79).  Furthermore,  we  have 
the  theorem  below. 


Theorem  Al.  (5j*'\  5®)  is  an  isolated  global  minimizer  of  G\{si,  S2\  0). 


Proof  of  theorem  Al.  In  the  proof,  A  is  fixed  at  A  =  0.  We  consider  the  free  energy  of  pure 
nematic  (non-dipolar)  polymers  for  probability  density  of  the  form  (78): 


G2(5i,  52)  = 


pirn)  In  p{m)dm  —  -{m  ^ 


(82) 


In  [11],  it  was  concluded  that  for  b  >  7.5  only  the  family  of  prolate  equilibrium  states  is 
stable  with  respect  to  the  Maier-Saupe  interaction  and  is  the  global  minimizer  of  the  free 
energy  of  the  pure  nematic  (non-dipolar)  polymers.  The  isotropic  solution  and  the  family  of 
oblate  solutions  are  unstable  for  b  >  7.5.  In  other  words,  the  probability  density  (78)  with 
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^2°^)  is  a  global  minimizer  of  the  right-hand  side  of  (82).  Consequently,  ^2°^)  is  a 
global  minimizer  of  G2(ii,  ^2)-  It  should  be  pointed  out  that  in  the  discussion  of  stability  we 
must  always  specify  what  are  the  interactions.  For  example,  a  stable  state  with  respect  to  only 
the  Maier-Saupe  interaction  may  be  unstable  with  respect  to  the  Maier-Saupe  interaction  plus 
the  dipole-dipole  interaction.  Using  form  (78),  we  can  write  G2{s\,  S2)  as 

G2{s\,S2)  =  b{si{m\)  +  S2{ml)  H-  (1  -  -  S2){ml))  -  In  Z  -  +  {mlf  -I-  {mlf  ) 


=  -  Ji)^  +  ({ml)  -  S2)^  +  ({ml)  -  (1  -  5i  -  ^2))^]  +  01(^1,52;  0) 


(83) 


which  implies 


Gi(ii,  S2',  0)  ^  G2(^l,  i2), 


Gi(ii,  ^2;  0)  =  G2(^i,  52)  if  (ji,  i2)  is  a  solution  of  (79).  (84) 

Since  (sf'\  i®)  is  a  solution  of  (79)  and  a  global  minimizer  of  G2(si,  ^2),  it  follows  that 
is  a  global  minimizer  of  Gi(si,S2',  0).  As  we  pointed  out  earlier,  all  global 
minimizers  of  Gi(si,  S2\  0)  are  solutions  of  equation  (79).  Property  (84)  also  tells  us  that 
for  b  >  7.5  neither  the  isotropic  state  nor  the  oblate  equilibrium  state  is  a  global  minimizer 
of  Gi  (i'l ,  S2',  0)  because  neither  of  them  is  a  global  minimizer  of  G2(s\ ,  82)-  Suppose  (si ,  S2) 
corresponds  to  the  isotropic  state  or  an  oblate  equilibrium  state.  Since  (i"i,  S2)  is  a  solution 
of  (79),  we  have  Gi(si,  S2,  0)  =  G2(il,  S2).  For  b  >  7.5,  a  prolate  equilibrium  state  is  a 
global  minimizer  of  the  free  energy  while  neither  the  isotropic  state  nor  an  oblate  equilibrium 
state  is.  Thus,  we  must  have  G2 (si,  52)  >  G2(5i*'\  52°^).  It  follows  that  G 1  (si ,  i2 >  0)  > 
Gi(s{°\  S2°\  0).  That  is,  (s”i,s'2)  is  not  a  global  minimizer  of  Gi(si,S2;  0).  Therefore,  a 
global  minimizer  of  Gi(si,  S2;  0)  must  come  from  the  family  of  rotation-equivalent  prolate 
equilibrium  states.  In  the  framework  of  probability  density  (78),  the  principal  axes  of  the 
second  moment  are  restricted  to  be  parallel  to  the  axes  of  the  coordinate  system.  With  this 
restriction,  the  family  of  rotation-equivalent  prolate  equilibrium  states  consists  of  only  three 
solutions:  (sj°\  S2°^),  (s®,Sj°^)  and  (s[°\  s[  ').  Therefore,  (s[°\  s®)  is  an  isolated  global 
minimizer  of  Gi(si,  S2;  0).  This  completes  the  proof  of  theorem  Al. 

Let  R(d)  denote  the  rectangle  [s®  —  d,  s®  -i-fif]  x  [s®  —  d,  s®  +d]  and  let  dR(d)  denote 
the  boundary  of  the  rectangle.  Theorem  Al  tells  us  that  there  exists  do  >  0  such  that 

Gi(ii,  ^2;  0)  >  Gi(^f ,  if;  0)  for  (^1,^2)  e  R(do)\(sf\ if).  (85) 

Let  S  =  min(i/o,  where  i^°^  >  0  is  the  order  parameter  of  the  prolate  solution.  Let 

D(X)=  min  Gi(ii,i2;  A.)  -  Gi(if  ,if ;;.).  (86) 

(si,s2)edR(S) 


Property  (85)  implies  that  D(0)  >  0.  Differentiating  (81)  with  respect  to  X  yields 


—  Gi(ii,  i2;  X) 
aX 


=  |{m3)|  <  1. 


We  select  Xi  =  D(0)/4  >  0.  Since  5jGi(ii,i2;  A)  is  uniformly  bounded  by  1,  we  have 
D(A.i)  ^  D(0)/2  >  0,  which  implies  that  function  Gi(ii,  i2;  Ai)  must  attain  a  minimum 
inside  the  rectangle  R(S).  Recall  that  all  critical  points  of  Gi(ii,i2;  X)  are  solutions  of 
equation  (79).  Thus,  for  Ai  =  D(0)/4  >  0,  equation  (79)  has  a  solution  (if,  sf )  e  R(S). 
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Once  we  know  ^2'^),  equation  (77)  is  satisfied  by  selecting  ai  —  X\/ {m^)  >  0.  In  this 
way,  we  constructed  an  equilibrium  solution.  From  the  selection  of  &  given  above,  we  obtain 


ic(0) 

2^ 


<  i 


(1) 


<  I 


1.(0) 

6^ 


1  +  1^(0) 

3  -r  2^ 


<  S. 


(1) 


2.(0) 

3^ 


.(1) 


<  i  +  55W 

^  3  T  6^ 

=  1  -  - 


(87) 


.(1) 


<  I 


where  >  0  is  the  order  parameter  of  the  prolate  equilibrium  solution.  Equation  (87) 
shows  that  '  is  still  the  only  eigenvalue  (of  the  second  moment)  that  is  above  1  /3.  Both 
and  53^^  are  below  1/3.  If  the  equilibrium  solution  were  to  be  axisymmetric,  then  the  only 
possible  choice  for  the  axis  of  symmetry  would  be  the  y-axis.  As  we  pointed  out  earlier  in  the 
discussion  of  equation  (77),  for  1  >  0  we  have  (m3)  >  0  which  implies  that  the  first  moment 
(polarity  vector)  is  parallel  to  the  z-axis.  Therefore,  it  is  clear  that  the  equilibrium  solution  we 
constructed  above  is  non-axisymmetric. 
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